6 research outputs found
Hiking in the energy landscape in sequence space: a bumpy road to good folders
With the help of a simple 20 letters, lattice model of heteropolymers, we
investigate the energy landscape in the space of designed good-folder
sequences. Low-energy sequences form clusters, interconnected via neutral
networks, in the space of sequences. Residues which play a key role in the
foldability of the chain and in the stability of the native state are highly
conserved, even among the chains belonging to different clusters. If, according
to the interaction matrix, some strong attractive interactions are almost
degenerate (i.e. they can be realized by more than one type of aminoacid
contacts) sequence clusters group into a few super-clusters. Sequences
belonging to different super-clusters are dissimilar, displaying very small
() similarity, and residues in key-sites are, as a rule, not
conserved. Similar behavior is observed in the analysis of real protein
sequences.Comment: 17 pages 5 figures Corrected typos added auxiliary informatio
Simple models of protein folding and of non--conventional drug design
While all the information required for the folding of a protein is contained
in its amino acid sequence, one has not yet learned how to extract this
information to predict the three--dimensional, biologically active, native
conformation of a protein whose sequence is known. Using insight obtained from
simple model simulations of the folding of proteins, in particular of the fact
that this phenomenon is essentially controlled by conserved (native) contacts
among (few) strongly interacting ("hot"), as a rule hydrophobic, amino acids,
which also stabilize local elementary structures (LES, hidden, incipient
secondary structures like --helices and --sheets) formed early
in the folding process and leading to the postcritical folding nucleus (i.e.,
the minimum set of native contacts which bring the system pass beyond the
highest free--energy barrier found in the whole folding process) it is possible
to work out a succesful strategy for reading the native structure of designed
proteins from the knowledge of only their amino acid sequence and of the
contact energies among the amino acids. Because LES have undergone millions of
years of evolution to selectively dock to their complementary structures, small
peptides made out of the same amino acids as the LES are expected to
selectively attach to the newly expressed (unfolded) protein and inhibit its
folding, or to the native (fluctuating) native conformation and denaturate it.
These peptides, or their mimetic molecules, can thus be used as effective
non--conventional drugs to those already existing (and directed at neutralizing
the active site of enzymes), displaying the advantage of not suffering from the
uprise of resistance
Deriving amino acid contact potentials from their frequencies of occurence in proteins: a lattice model study
The possibility of deriving the contact potentials between amino acids from
their frequencies of occurence in proteins is discussed in evolutionary terms.
This approach allows the use of traditional thermodynamics to describe such
frequencies and, consequently, to develop a strategy to include in the
calculations correlations due to the spatial proximity of the amino acids and
to their overall tendency of being conserved in proteins. Making use of a
lattice model to describe protein chains and defining a "true" potential, we
test these strategies by selecting a database of folding model sequences,
deriving the contact potentials from such sequences and comparing them with the
"true" potential. Taking into account correlations allows for a markedly better
prediction of the interaction potentials
Multiple Routes and Milestones in the Folding of HIV–1 Protease Monomer
Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease
Energy Profile of the Space of Model Protein Sequences
A numerical study of the energy landscape of the space of model protein sequences is carried out. As a consequence of the heterogeneity of the contact energies among amino acids, the energy landscape displays a very rough profile, a behaviour typical of frustrated systems. This gives raise to a hierarchical clustering of low-energy sequences and can have evolutionary consequences